Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 February 2008 (MN WF$t style file v2.2) 



Large-Scale Non-Gaussianities in the 21 cm Background 
Anisotropics From the Era of Reionization 

Asantha Cooray 

Theoretical Astrophysics, MS 130-33, California Institute of Technology, Pasadena, CA 91125. 
Dept. of Physics and Astronomy, 4129 Frederick Reines Hall, University of California, Irvine, CA 92697. 
E-mail: asante@caltech.edu 

2 February 2008 

ABSTRACT 

The brightness temperature fluctuations in the 21 cm background related to neutral 
Hydrogen distribution provide a probe of physics related to the era of reionization 
when the intergalactic medium changed from a completely neutral to a partially ion- 
ized one. We formulate statistics of 21 cm brightness temperature anisotropies in 
terms of the angular power spectrum, the bispectrum, and the trispectrum. Using the 
trispectrum, we estimate the covariance related to the power spectrum measurements 
and show that correlations resulting from non-Gaussianities are below a percent, at 
most. While all-sky observations of the 21 cm background at arcminute-scale resolu- 
tion can be used to measure the bispectrum with a cumulative signal-to-noise ratio 
of order a few ten, in the presence of foregrounds and instrumental noise related to 
first-generation interferometers, the measurement is unlikely to be feasible. For most 
purposes, non-Gaussianities in 21 cm fluctuations can be ignored and the distribution 
can be described with Gaussian statistics. Since 21 cm fluctuations are significantly 
contaminated by foregrounds, such as galactic synchrotron or low-frequency radio 
point sources, the lack of a significant non-Gaussianity in the signal suggests that any 
significant detection of a non-Gaussianity can be due to foregrounds. Similarly, in ad- 
dition to frequency information that is now proposed to separate 21 cm fluctuations 
from foregrounds, if the non-Gaussian structure of foregrounds is a priori known, one 
can potentially use this additional information to further reduce the confusion. 

Key words: cosmology: theory — diffuse radiation — large scale structure — inter- 
galactic medium 



1 INTRODUCTION 



The 21 cm radio background related to spin-flip transition 
is a well-known probe of the neutral Hydrogen density con- 
tent prior to and during the era of reionization (Scott & 
Rees 1990; Madau et al. 1997; Zaldarriaga, Furlanetto & 
Hernquist 2004; Morales & Hewitt 2004). This background 
has the advantage that by selecting the observed frequency 
of the redshifted line emission, one can directly view the 
neutral Hydrogen distribution of the Universe at a given 
redshift. With maps of the neutral Hydrogen content, as 
a function of the redshift, one can obtain statistics related 
to the reionization process, especially during the transition 
from a fully neutral medium to an ionized one (Barkana & 
Loeb 2001). Similarly, prior to the epoch of reionization 
when the Universe is fully neutral, 21 cm fluctuations pro- 



vide a probe of the primordial density perturbations (Loeb 
& Zaldarriaga 2004; Bharadwaj & Ali 2004). 

In addition to 21 cm data, information related to the 
era of reionization can be obtained with large angular scale 
cosmic microwave background (CMB) polarization observa- 
tions. With CMB anisotropies, one is studying the pres- 
ence of electrons through various Compton-scattering re- 
lated processes where modifications are introduced to the 
temperature and polarization fluctuations. Unlike 21 cm ob- 
servations, unfortunately, CMB data only allow a measure- 
ment of the integrated column density of electrons (Zaldar- 
riaga 1997) with no detailed information on the exact reion- 
ization history of the universe (Hu & Holder 2003; Kapling- 
hat et al 2003). As has been studied recently, Lyman-a 
optical depth from Gunn-Peterson troughs (Gunn & Peter- 
son 1965) towards z ~ 6 quasars in the Sloan Digital Sky 
Survey (Fan et al 2002) provides information related to 
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the end of reionization process when the Universe is almost 
fully ionized. While recent measurements from the WMAP 
mission (Bennett et al 2003) suggest an optical depth to 
electron scattering of 0.17 ± 0.04 (Kogut et al 2003), it 
can be explained with a variety of reionization models that 
also take into account the ionizing fraction estimates from 
z ~ 6 quasars (Cen 2003; Ciardi & Madau 2003; Chen et 
al 2003; Gnedin 2004; Wyithe & Loeb 2003; Haiman & 
Holder 2003). 

With 21 cm observations, one can potentially recon- 
struct the exact reionization history and determine, espe- 
cially, if the reionization process was either instantaneous 
or lengthy and complex as implied by current CMB and 
Sloan data. Moreover, since observable effects in CMB data 
depend only on the ionized content while 21 cm fluctua- 
tions come from inhomogeneities in the neutral distribution, 
the combination is expected to allow additional information 
related to the era of reionization (Cooray 2004; Cooray 
2004). Detailed analytical and numerical studies related to 
the physics associated with the 21 cm background are now 
motivated by progress in the experimental side involving 
plans for a variety of low-frequency radio experiments such 
as PAST 1 (Pen 2004), Low Frequency Array (LOFAR) 2 , 
and Square Kilometer Array (SKA) 3 . 

While the scientific motivation to study 21 cm fluctua- 
tions during the era of reionization is strong, there are var- 
ious challenges to overcome. These include the removal of 
foreground radio emission, such as due to synchrotron emis- 
sion from the Milky Way (Shaver et al 1999), low frequency 
radio point sources (Di Matteo et al 2002) and free-free 
emission from free electrons in the intergalactic medium (Oh 
1999; Cooray & Furlanetto 2004). Contrary to initial sug- 
gestions that 21 cm background studies may be impossible 
due to the large number of contaminating sources (Oh & 
Mack 2003), with brightness fluctuations much higher than 
expected for neutral Hydrogen at high redshifts, recent sug- 
gestions based on multifrequency analysis indicate that fore- 
ground contamination can be reduced given the smoothness 
of these confusions in frequency space (Shaver et al 1999; 
Di Matteo et al 2002; Zaldarriaga, Furlanetto & Hernquist 
2004; Santos, Cooray & Knox 2004). 

In current studies, the statistics related to 21 cm fluc- 
tuations have been mostly restricted to the angular power 
spectrum of anisotropies. This is primarily motivated by the 
expected Gaussian nature of these fluctuations. It is useful to 
consider if there are non-Gaussian signals in the 21 cm back- 
ground anisotropies and, if any exists, whether these signals 
can be measured with planned low-frequency instruments. 
Here, we discuss the role of non-Gaussianities in the 21 cm 
background. For this purpose, we describe spatial fluctua- 
tions in the 21 cm background during the era of reionization 
as due to a combination of spatial fluctuations in the under- 
lying neutral Hydrogen density, peculiar velocity of neutral 
gas, and the fraction of ionized gas relative to the total gas 
density. The non-Gaussian signals are generated by the non- 

1 http://astrophysics.phys.cmu.edu/ jbp 

2 http://www.lofar.org 

3 http://www.skatelescope.org 



linear coupling of fluctuations in either the neutral density 
and peculiar velocity or neutral density and ionized-fraction 
fields. 

We model the trispectrum, the Fourier analog of the 
four-point correlation function, of 21 cm anisotropies and ex- 
plore how these non-Gaussianities affect the power spectrum 
measurement in terms of its full covariance. The trispectrum 
produces an additional source of noise beyond the stan- 
dard cosmic variance under Gaussian statistics and corre- 
lates power spectrum measurements between different mul- 
tipoles (Scoccimarro, Zaldarriaga & Hui 1999; Cooray & 
Hu 2001b). As we find, in the case of the 21 cm background, 
these correlations are at the level of at most a percent sug- 
gesting that non-Gaussianities are not likely to be a ma- 
jor source of complication when making a power spectrum 
measurement. We also estimate the bispectrum, the three- 
point correlation in Fourier space, of 21 cm anisotropies. 
While all-sky observations with no significant instrumental 
noise can be used to measure the bispectrum with signal- 
to-noise ratios of a few ten, for more practical observations 
with the first-generation of low-frequency radio interferom- 
eters, the signal-to-noise ratio is below a unity. In general, 
21 cm background anisotropies can be described with Gaus- 
sian statistics and measurements can be focused primarily 
on the power spectrum. In our calculations, we mainly focus 
on the large physical, and thus large angular, scale fluctu- 
ations of the 21 cm background. Non-linear effects related 
to the neutral gas distribution and velocities can lead to 
additional non-Gaussian signals, but such effects are best 
captured in numerical simulations of the reionization pro- 
cess. These include the effects associated with finite size of 
reionized regions (c.f., Bharadwaj & Pandey 2004), and we 
leave these and other issues to be discussed in future works. 
For first-generation observations, such effects are unlikely to 
be important given restrictions on the angular resolution of 
observations. The lack of a significant non-Gaussianity at 
large angular scales, at least, has the advantage that if any 
significant non-Gaussianity is detected, it can be ascribed to 
foreground sources. Thus, in addition to smooth behavior in 
frequency space, the non-Gaussian structure of foregrounds, 
such as those related to Galactic synchrotron, can be used 
to further reduce the confusion. 

The paper is organized as follows: in we discuss the 
21 cm signal anisotropy, the trispectrum in terms of the 
power spectrum covariance, and the bispectrum. While we 
only consider the non-Gaussian structure of the 21 cm back- 
ground, the lack of a significant non-Gaussianity in the sig- 
nal can be potentially used to separate nearby galactic fore- 
grounds which are expected to be highly non-Gaussian. We 
discuss these and other possibilities in JJ] Throughout the 
paper, we make use of the WMAP-favored ACDM cosmo- 
logical model (Spergel et al 2003). 



2 THE 21CM SIGNAL 

The 21 cm surface brightness temperature observations are 
primarily observed as a change in the intensity of CMB due 
to line emission and absorption. The change in brightness 
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Figure 1. Top: The power spectrum for the 21cm signal at z=10.0 
and assuming a bandwidth of Su = 1MHz. Units are in mK 2 . 
The calculation makes use of a reionization model from Santos, 
Cooray & Knox (2004), though similar results are obtained for 
slightly varying models from Santos et al. (2003). Here, we show 
the total contribution to the power spectrum (solid line), as well 
as the power spectrum related to neutral density (dotted line) 
and peculiar velocity (dashed line) fluctuations alone. We also 
show the Gaussian cosmic variance (long-dashed line) and an ad- 
ditional contribution to the power spectrum variance related to 
non-Gaussian trispectrum (dot-dashed line; see §2.1 for details). 
The solid line labeled 'noise' shows an estimated instrumental 
noise curve related to 21 cm observations with the SKA following 
Zaldarriaga, Furlanetto & Hernquist (2004). Bottom: The ratio of 
various contributions to the 21 cm anisotropy power spectrum. In 
the top curves, we show the ratio of total contribution to that re- 
lated to density fluctuations alone, while the bottom set of curves 
are the ratio of peculiar velocity contribution to that of the total 
power spectrum. In each of these two sets, the four curves, from 
top to bottom, are for observational bandwidths of 0.1, 0.5, 1.0 
and 2.5 MHz, respectively. 



temperature, T2i(v), compared to the CMB at an observed 
frequency v is then 



Ts — Tcmb 

1+2 



(1) 



where Ts is the spin temperature of the IGM, z is the red- 
shift corresponding to the frequency of observation (1 + z = 



vsi/v, with u 2 i = 1420 MHz) and T C mb = 2.73(1 + z)K is 
the CMB temperature at redshift z. 

In order to simplify the calculation related to 21 cm fluc- 
tuations, we assume that T s ^> Tcmb, in a scenario where 
the gas density field in the intergalactic medium is heated as 
part of the reionization and structure formation processes. 
These include shocks associated with the formation of virial- 
ized dark matter halos where first stars are eventually found, 
and whose UV photons subsequently reionize the universe, 
and X-rays from first supernovae (Venkatesan, Giroux & 
Shull 2001; Chen & Miralda-Escude 2004). At sufficiently 
high z, prior to significant reionization and heating, we ex- 
pect T a < Tcmb, and the HI density content is seen as an 
absorption in 21 cm with respect to the CMB (Madau et 
al. 1997). The statistics of 21 cm fluctuations related to 
such an early epoch can be used as a probe of primordial 
fluctuations (Loeb & Zaldarriaga 2004; Bharadwaj & Ali 
2004). 

The optical depth, r, of this patch in the hyperfine tran- 
sition (Field 1959) is given in the limit of fesT s >> hv2i 
by 



3c hAio nm 
Vokvl x T s H{z) 



(2) 
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where A10 is the spontaneous emission coefficient for the 
transition (2.85 x 10 -15 s _1 ), nm is the neutral Hydrogen 
density, and v is the peculiar velocity of the neutral gas dis- 
tribution when r is the comoving radial distance and H(z) 
is the expansion rate at a redshift of z. The neutral density 
be expressed as nm = xnn g (l + S g ), when n g is the mean 
number density of cosmic baryons, with a spatially varying 
overdensity S g and xh is the fraction of neutral Hydrogen 
(xh = 1 — x e where x e is the fraction of free electrons). 
We refer the reader to Zaldarriaga, Furlanetto & Hernquist 
(2004) and Santos, Cooray & Knox (2004) for further de- 
tails. Note, however, that prior calculations related to the 
21 cm anisotropies ignored both peculiar velocities and ion- 
izing fraction fluctuations (e.g., Zaldarriaga, Furlanetto & 
Hernquist 2004), or peculiar velocities (Santos, Cooray & 
Knox 2004). 

With T s » Tcmb, the 21 cm signal will be observed 
in emission with respect to the background CMB. Writing 
the brightness temperature as: 



r 2 i(z) = c [i + Sg- 

where 
cs; 23 



z dv 



(1 — X e — X e 5 x ) , 



(3) 



0.02 



H(z) dr 



and 



8 X is the perturbation in the ionization fraction (S x = 
x e -x e y Here and above, the spatial dependent quantities 
are evaluated at the position x = rn and for simplicity in 
notation we have dropped this dependence. 
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The measured brightness temperature corresponds to a 
convolution of the intrinsic brightness with some response 
function Wv{r) that characterizes the bandwidth of the ex- 
periment: 



T(n, vo) 



drW„ (r)T 21 {nr,r), 



(•») 



where n is the direction of observation and r corresponds to 
the observed frequency v. We can now determine the angular 
power spectrum of T(n, i/q) which is related to the 3-d power 
spectrum of T2i(n, r), defined by 

(f 2 i(k,u 1 )f 2 i{k',u 2 )) = (2irfS D (k + k')P 21 (k,» 1 ,v 2 ), (6) 



where T 2 i(n,r) = T 21 (r, v) = / ^fai (k, i/)e lk ' r . 

Making use of the spherical harmonic moment of the 21 
cm fluctuations, at a frequency vq, 



ai m {v )= / dnY^(n)T(n, u ) , 



(7) 



we define the angular power spectrum as (Kaiser 1992) 

(Oi m (l/i)ofm(l^)) = Ci(h>i,V 2 ) . (8) 

To calculate this power spectrum, first, we write Fourier 
space fluctuations in the brightness temperature as 

T 21 (k,z) = c[8 g (k)(l-x e )+^ 2 8 v (k)(l-x e )~x e S x (k) 
-■>: I 4-^^(k')5 g (k-k') 



(2tt) 3 



-X e fl 



tfk' 

(2tt) s 



^(k')^(k-k')] (9) 



where 5 V is the velocity divergence, \i = n • k, and // = 
n • (k — k'). Here, we concentrate primarily on the large 
scale velocity fluctuations and at these scales, we can use 
linear theory to write S v {k) = /(Sl m )5 lln (k) where /(f2 OT ) = 
91nG/91naw Q^ 6 + 1/70[1 - l/2fi m (l + Q m )], where G is 
the growth factor related to evolution of density perturba- 
tions with the approximation valid in a flat cosmology with 
a cosmological constant. 

To calculate the power spectrum of 21 cm fluctuations, 
we ignore terms which are second order and higher and 
involving the convolution between density and ionization- 
fraction and velocity and ionization-fraction fluctuations. 
We make use of the Rayleigh expansion for the plane wave 
given by 



e*- = 4nJ2i'ji(kr)Yr(AWr(i) , 



(10) 



and separately consider fluctuations in the velocity and den- 
sity fields separately. To simplify the calculation associated 
with the velocity fluctuations, we also use 



n-k = ^^r(fl)C(k). 



(11) 



After straightforward, but tedious, algebra, we can 
write the power spectrum as 

Ci(ui, v 2 ) = — / k 2 dk 

7T / 



x [p aa (k, v x ,v 2 )I^ (k)Ip(k) + P bb (k, ui,u 2 )Jp (k)Jp (fc) 
+ P ah {k,v^ 2 )[I»\k)J^(k) + J^(k)I^{k)]} (12) 



where 



drWv(r)ji(kr) 



drW v {r)j'i \kr) 



(13) 



with the spherical Bessel function given by ji(kr). In deriv- 
ing Eq. 1121 we have ignored terms which are second order 
in 8. In Eg. 1131 the double derivative of the spherical Bessel 
function is with respect to its argument. 

The corresponding three-dimensional power spectra are 

P aa (k,z) = c [{1- x e ) 2 P gg {k,z) +x 2 e Ps I s !C {k,z) 



(14) 



Pbb(k,z) 

Pab(k,z) 



2-Pgix (k, z)x e (l — x e )] , 
2 [(1 - x e ) 2 Pg v s v (k, «)] 



[(1 - x e ) 2 P gSv (k, z) - P Sv s x (k, z)x e (l - x e jj 



The neutral Hydrogen density power spectrum is repre- 
sented by P gg (k,z), Ps x s x (k,z) is the power spectrum re- 
lated to perturbations in the ionized fraction, Pgs^. (k, z) 
is the cross-correlation power between these two, and 
Ps v s v (k, z) is the velocity divergence spectrum. Using linear 
theory evolution of gravitational perturbations, the latter 
can be written as P$ v $ v (k, z) fa f 2 (fl m )P hn (k,z) in terms 
of the linear density spectrum. The cross spectra between 
velocity and neutral density and ionized fraction are given 
as P g s v and Ps v s x {k, z), respectively. 

We assume the neutral Hydrogen density field is an 
unbiased tracer of the dark matter density field and write 
Pgg(k,z) = Pss(k,z) and, when illustrating our results, use 
the linear theory spectrum to describe dark matter around 
the era of reionization. More complicated models that cap- 
ture the non-linear behavior of the gas density field can be 
obtained through approaches such as the halo models of 
large scale structure (Cooray & Sheth 2002), but extended 
to higher redshifts. Note that the correlations in the ionized 
fraction depend on the reionization history and we make use 
of the patchy reionization model in Santos et al (2003) to 
write 

PsA {k, z) = b 2 (z)P ss (k, z)e~ k2R2 (15) 

where b(z) is a mean bias that captures the halo bias (Mo 
et al. 1997) weighted by the different halo properties, and 
R the mean radius of the HII patches. The time-dependence 

of R is modeled as R = ( ^jg ) 1,/3 R P where R p is the (co- 
moving) size of the fundamental patch taken to be ~ 100 
Kpc. 

In addition to the power spectra of the gas density, ve- 
locity divergence, and ionization fraction, when presenting 
numerical results, we assume cross-correlation spectra be- 
tween each of these three fields can be described with perfect 
correlations such that for example, P xy — y/ P xx P vy where 
x and y are any of the two fields of interest. 

In Fig. 1, we show the power spectrum of 21 cm fluctu- 
ations at a redshift of 10 with a bandwidth for observations 
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Figure 2. The correlation coefficients related to the non- 
Gaussian covariance of 21 cm power spectrum measurements at 
z = 10 with a bandwidth for o bservations of 1 MHz. Here, we 
show Covj^ / y/ CovY^i CovJ^ as a function of I2 when l\ = 10, 
100, and 1000. Note that these are correlations only in the non- 
Gaussian part of the contribution to the covariance. The total 
correlations, when considering the measurement of the 21 cm an- 
gular power spectrum, include the part due to Gaussian variance. 
Since Gaussian variance is roughly two orders of magnitude or 
more higher than that of non-Gaussianitics, as shown in Fig. 1, 
these correlations are at most a percent and can be ignored for 
all practical purposes. 



of 1 MHz. Here, we make use of a reionization model from 
Santos, Cooray & Knox (2004), but similar results, with 
amplitudes varying as shown in Fig. 2 of Santos, Cooray 
& Knox (2004), are found for other models. The reioniza- 
tion model provides information on x e and ionizing source 
bias. In Fig. 1, we separate the contribution to the 21 cm 
anisotropy power spectrum to those related to the neutral 
density /ionizing fraction fields (P aa ) and peculiar veloci- 
ties (Pbb)- In addition to these two, the cross spectra be- 
tween these quantities {P a b) also contribute to the total 21 
cm anisotropy angular power spectrum and this leads to a 
roughly a factor of 2 increase in power at large angular scales 
when compared to density fluctuations alone. 

In the bottom panel of Fig. 1, we show how this fraction 
varies as a function of the observational bandwidth. In ad- 
dition to the ratio of total to density contributions, we also 
show the fractional contribution related to peculiar velocity 
fluctuations alone. For all bandwidths, this contribution is at 
the level around 0.3; one cannot simply ignore the presence 
of peculiar velocities based on assumptions related to a large 
bandwidth for observations (c.f., Zaldarriaga, Furlanetto & 
Hernquist 2004) . Note that the importance of peculiar veloc- 
ities has already been discussed in Bharadwaj & Ali (2004) 
in the context of 21 cm fluctuations prior to reionization 
in the form of an absorption signal related to background 
CMB. 



2.1 Power spectrum covariance 

In prior works related to the 21 cm anisotropies, the fluc- 
tuations were assumed to be Gaussian so that the power 
spectrum captures all statistical information and its vari- 
ance is simply given by Gaussian statistics (e.g., Zaldarriaga, 
Furlanetto & Hernquist 2004; Santos, Cooray & Knox 2004). 
The second order corrections to the 21 cm brightness tem- 
perature, in an era of partial reionization, however, leads to 
a non-Gaussian contribution to the 21 cm brightness tem- 
perature anisotropies. These non-Gaussianities, then, add an 
extra covariance component to the power spectrum measure- 
ment and a reliable estimate of the signal-to-noise ratio of 
the power spectrum measurement should account for these, 
and other potentially interesting, non-Gaussian effects. 
Writing the measured power spectrum as 



C, 



21+1 ^ 

m— — l 



(aim(Hm) ■ 



(16) 



we note that the full covariance of this measured power spec- 
trum is 

Cov ili2 = (17) 

21 + 121 + 1 ^ ai ^ m ^ a h rn i ai 2 m 2 a hm 2 ) ~ Cl x Cl 2 , 

mi m 2 

and separate this to part which is Gaussian and a part which 
is non- Gaussian such that 



Cov, li2 = Cov^ 2 u + CovE? 2 



(18) 



Cov 



Cov, 



Cfi<5ii! 



2h + l 



2ll + 1 212 + 1 ^ ' (o-hmx dl imi a l2m 2 a l 2 m 2 ) c 



where (..) c denotes the connected part of the four-point cor- 
relator. To calculate this non-Gaussian part, we now make 
use of the second order correction and, here, only focus on 
the part involving the convolution between neutral Hydro- 
gen density and ionization fraction fluctuations. In terms of 
multipole moments, this second order term can be written 
as 



a lr. 



-(4nf 



(2 



|ki_ f 



)3 y, y, 

h m l l 2 m 2 



H) 



h + h 



x J drcx e ji 1 (k 1 r)ji 2 (k 2 r)S x (k 1 )S g (k2) (19) 

x ^(kO^ 2 *^) / dhYr(n)Y^(n)Y^(h). 
Using this term twice, we simplify to obtain 



imi Q>1 imiOl 2 m 2 0* 2m2 >c 



(20) 



'i\dki / k\dki I k^dks.c'xi 



x 2 e P XeXe (fci)P XeXe (k 2 )P gg (k 3 ) 
+ (l-x e ) 2 P 99 (fc 1 )P 99 (fc 2 )P^ e (fc 3 ) 
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X ^ Il 2 {k2)Il 1 (kl)Kl 1 ,l(kl,k3)Kl 2 ,l(k2,k 3 ) 

rn\m,2lm 



X 

where 

K h A 2 (k 1 ,k 2 ) 



1 



drW„(r)ji 1 (k 1 r)jt 2 (k 2 r) . 



(21) 



Including permutations, the non-Gaussian part of the power 
spectrum covariance is then 

CovfJ? 2 = (22) 



/ kldk^xl 



x e P XcX .(k 1 )P XcXc {k2)P gg {h A ) 



+ (l-x e rP 99 (fci)P 99 (fc 2 )P^ e (fc 3 ) 

x Ii 2 (k2)Ii 1 (k 1 )Ki 1:t (k 1 ,k 3 )Ki 2t i(k 2 ,k 3 ) , 

and an additional permutation with the replacement of 
h — ► £2. Here, in simplifying, we have made use of the 
orthonormality relation between Wigner-3j symbols. Note 
that we have ignored the cross power spectra between den- 
sity and ionization fraction here since these are likely to 
be subdominant relative to the density and ionization frac- 
tion power spectra alone. Furthermore, we have ignored the 
non- Gaussian signals associated with the coupling between 
density fluctuations and peculiar velocities. These are sub- 
dominant due to mismatch in projection between density 
and velocity fluctuations, captured by projection integrals 
of Ii with a ji, and J; with a j" . 

In Fig. 1, we show the variance related to Gaussian 
cosmic variance term and the non-Gaussian contribution 
separately in the form of an error on the measured power 
spectrum and unbinned in multipole space. These errors are 
equivalent to that of \/Covu but weighted by a factor of 
1(1 + l)/27r. As shown, the non-Gaussian correction to the 
unbinned power spectrum measurement error is two orders 
of magnitude below when compared to that due to Gaussian 
errors. 

In Fig. 2, we show the correlation coefficient related 
to the no n-Gaussian cov ariance where we plot Corr — 
Covf ± f 2 1 Cov^Y Cov^ 2 as a function of l 2 for several val- 
ues of h. While these correlations are large, note that what 
matters for the power spectrum measurement is the to- 
tal correlations, or rather, Cov^ / yj Cov;^ CoV( 2 ; 2 . Since 
Cov hh = Cov°?"+Cav£° ~ Cov^j", these correlations are 
below a percent at most (when I ~ few 10 3 ) and can be ig- 
nored for all reasonable purposes involving the measurement 
of the 21 cm power spectrum. 



2.2 Bispectrum 

While corrections to the power spectrum covariance from 
non-Gaussianities are negligible, it is still useful to con- 




Figure 3. Cumulative signal-to-noise ratio for the measurement 
of the 21 cm bispectrum as a function of l 3 , for observations 
at z = 10 with a bandwidth of 1 MHz. From top-to-bottom, the 
three curves for an all-sky experiment with no instrumental noise, 
a 4000 deg. 2 survey with instrumental noise as show in Fig. 1, and 
a survey of 400 deg. 2 and with an increase in noise by a factor 
of 2 to account for unsubtracted foregrounds following Santos, 
Cooray & Knox (2004). While the signal-to- noise ratio is at the 
level of a few tens for the perfect experiment, for more realistic 
scenarios, the signal-to-noise ratio is below a unity. 



sider if non-Gaussianities 
measured with upcoming 
non-Gaussian signal, here, 
der statistic and estimate 
ated with the bispectrum, 
the three-point correlation 
point correlation function 



in the 21 cm background can be 
interferometers. To calculate the 
we concentrate on the lowest or- 
the signal-to-noise ratio associ- 
or the Fourier space analogue of 
function. For reference, the three- 
is related to multipole moments 



B(n,m,l) = <T(n)T(m)T(i)> (23) 
= ^{ai^ai^ai^Y™ 1 (n)Y™ 2 (m)Y™ 3 (I) , 

where the sum is over (h, mi), (l 2 , m 2 ), (h, m 3 ). Statistical 
isotropy allows us to express the correlation in terms an m- 
independent function, 



(ai imi ai 



im 2 tti 3 m 3 



3) 



h h h 
mi rn 2 m 3 



(24) 



where the quantity Bi^is i s described as the angular aver- 
aged bispectrum (Cooray & Hu 2000) and we have made use 
of the orthonormality relation related to Wigner-3j symbols. 
As discussed in Cooray & Hu (2000), the angular bispec- 
trum, Bi 1 i 2 ; 3 , contains all the information available in the 
three-point correlation function and frequently used quan- 
tities such as the skewness and the collapsed three-point 
function can be written in terms of the bispectrum, such as 
a filtered sum of certain bispectrum configurations (Cooray 
& Hu 2001a). 

We consider the bispectrum formed by the 21 cm bright- 
ness fluctuations. This non-Gaussianity arises from the sec- 
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ond order terms in Eq. 3 (or in Fourier space, Eq. 9). Here, 
we again ignore the correction related to peculiar velocities, 
but instead focus on the second order term related to the 
coupling between neutral gas density fluctuations and the 
ionization fraction distribution. The bispectrum is simply 
then 



{ a ; L m! G^2 m 2 a h m 3 ) 



I kjdk, 



k 2 dk2c(l — x e )x 2 e Ii 2 {k\)Ii 3 (k2) 



x Ki 2 ,i 3 {kiM) [P gg (ki)P XeXe (k 2 ) + 
x J dnY^(h)Y t ^(h)Y^(n), 



(fci)JWfci)] 

(25) 



where, again, we have made use of the af m term from Eq. 19. 
The bispectrum is given by 



B h i 2 i 3 



E 



ii 

mi 



h 
m 2 



h 
m 3 



[(Vi 



at 



2 m 2 a l3m 3 



(2fi + 1)(2Z 2 + 1)(2Z 3 + 1) 
An 

x [bi 2 ,i 3 + Perm.] . 



h 




Here, 

b h,h 



= I fcidfei 



fc|dfac 3 (l — x e )x 2 



X [P gg (fa)P* e . e (fa) + Pgx e (kl)Pg Xe (fa)] 

xIi 2 (ki)Ii 3 (k2)Ki 2i i 3 (k 1 ,k 2 ) . 



(27) 



The permutations here account for the 3 additional terms 
with the replacement of fo — > h and l\ . 

The signal-to-noise ratio for the bispectrum is calcu- 
lated as 



/sky 



h>h>h 



^tot £*tot £*tot ' 
h l 2 h 



(28) 



Here, Cj ot represents all contributions to the power spec- 
trum of the ith field, 



^jtot ^ _|_ ^-tnoisc _|_ ^if'orcg 



^ , (29) 

where Cf° lsc is the instrumental noise contribution related 
to a 21 cm observations, and (7 ; forog ; s the confusing fore- 
ground contribution. The expression assumes all-sky obser- 
vations of the 21 cm background, but in the case of partial 
sky observations, the square of the signal-to-noise ratio is 
reduced by the fraction of sky covered by the data, / B ky 
We make use of the noise spectrum shown in Fig. 1 and 
include foregrounds in terms of the residual noise left over 
after multi-frequency cleaning following Santos, Cooray & 
Knox (2004). This procedure generally results in a factor of 
~ 2 increase in instrumental noise in the mutlipole range of 
general interest. 

In Fig. 3, we summarize our results related to the signal- 
to-noise ratio calculation. While the signal-to-noise ratio is 
at the level of a few tens for a no instrumental noise exper- 
iment that map out essentially the all sky, for partial sky- 
coverages and for added instrumental and foreground noise, 
the signal-to-noise ratio decreases below a unity suggesting 



that a reliable detection of the 21 cm anisotropy bispectrum 
may not be possible with first-generation low-frequency ra- 
dio interferometers. 



3 DISCUSSION & SUMMARY 

The brightness temperature fluctuations in the 21 cm back- 
ground related to neutral Hydrogen content provide a useful 
probe of physics related to the era of reionization during the 
transition from a completely neutral intergalactic medium to 
a partially ionized one. Here, we have presented a complete 
description of brightness temperature anisotropies related 
to 21 cm background in terms of fluctuations in the spa- 
tial distribution of density, peculiar velocity and ionization- 
fraction. Our calculations are more appropriate to the large 

+ J^erfrPT 9 ^' an< ^ ^ US ^ ar S e angular, scales since we have not at- 
tempted to model the small-scale non-linear behavior of the 
reionization process as well as effects associated with finite 
size of ionizing patches, among others. 

Our formalism allows us to calculate the angular power 
specteum, the bispectrum, the Fourier analog of the three- 
point correlation function, and the trispectrum, the Fourier 
analog of the four-point correlation function, of 21 cm fluc- 
tuations. In the case of the power spectrum, we have shown 
that peculiar velocities cannot easily be ignored since varia- 
tions in the large scale velocity field and its cross-correlation 
with the density distribution lead to an additional contribu- 
tion to the angular power spectrum of 21 cm anisotropies. 
For simple models related to the correlation between den- 
sity and velocity fluctuations, the correction is at the level 
of a factor of 2. As shown in Fig. 1, at large angular scales, 
this correction is independent of the bandwidth of 21 cm 
observations while with increasing bandwidth, the correc- 
tion reduces at arcminute angular scales corresponding to 
multipoles of a 10 3 and higher. 

We also extend previous discussions related to the 
power spectrum measurement by including an estimate of 
the non-Gaussian covariance resulting from parallelogram 
configurations of the trispectrum. In the case of the vari- 
ance, in terms of an unbinned error on the power spectrum, 
the non-Gaussian correction is at least two orders of magni- 
tude at degree angular scales with an increase to slightly less 
than an order of magnitude at multipoles of 10 4 (Fig. 1). 
In the case of the covariance, in terms of correlations be- 
tween power spectrum measurements at different multipoles, 
the correlations are at most at the level of a few percent, 
when I ~ 10 4 , suggesting that for most practical purposes 
non-Gaussianities can be safely ignored. While the effect of 
non-Gaussianities on the power spectrum measurement is 
insignificant, in terms of a direct measurement of the non- 
Gaussian signal, we have explore the possibility related to 
a measurement of the bispectrum. As shown in Fig. 3, the 
signal-to-noise ratio for the full bispectrum measurement is 
at most at the level of a few tens. For more practical ob- 
servations with the first generation experiments, this ratio 
drops below unity both due to added instrumental noise and 
the small sky-coverage. It is unlikely that the non-Gaussian 
signal will be directly measurable with 21 cm data alone. 
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It could be possible to pull out additional information, es- 
pecially related to the non-linear mode coupling second or- 
der terms, using cross-correlations between small scale CMB 
anisotropies (Cooray 2004b). 

Thus, for most practical purposes, the presence of non- 
Gaussianities can be ignored and the 21 cm fluctuations can 
be described with Gaussian statistics. Since 21 cm fluctua- 
tions are significantly contaminated with foregrounds, such 
as galactic synchrotron or low- frequency radio point sources, 
the lack of a significant non-Gaussianity in the signal sug- 
gests that any significant detection of non-Gaussianity can 
be easily ascribed to a foreground confusion. Similarly, in ad- 
dition to frequency information that is now proposed to sep- 
arate 21 cm fluctuations from foregrounds, if non-Gaussian 
structure of foregrounds is a priori known, one can poten- 
tially use this additional information to further reduce the 
foreground confusion. This requires a better understanding 
of the foreground distribution, especially involving the non- 
Gaussian nature of foregrounds. While in Santos, Cooray & 
Knox (2004), foregrounds were described only in terms of 
the power spectra, the formalism can be modified to include 
this additional information and we hope to return to this 
topic in a later discussion. 

While our calculations show that non-Gaussianities in 
the 21 cm background may not be significant, we may have 
underestimated the non-Gaussian nature by only concen- 
trating on the large physical, and thus large angular, scales 
of the signal. For example, additional non-Gaussianities are 
expected when the density field becomes non-linear. This, 
for redshifts z ~ 10, happens at scales below the typical res- 
olution of upcoming instruments and can be safely ignored. 
Similarly, the non-linear nature of the reionization process, 
especially the finite patch size of reionization patches, could 
lead to an additional non-Gaussian signal. This again de- 
pends on model dependent quantities and are better ad- 
dresses with numerical simulations rather than analytical 
calculations. 

Acknowledgments: Author acknowledges support from 
the Sherman Fairchild foundation and NASA NAG5-11985. 



Di Matteo, T., Perna, R., Abel, T., & Rees, M. 2002, ApJ, 564, 
576 

Fan X., Narayanan V. K., Strauss M. A., White R. L., Becker R. 

H., Pentericci L., Rix H. W. 2002, AJ, 123, 1247 
Field, G. B. 1959, ApJ, 129, 525 
Gnedin, N. Y. 2004, ApJ, 610, 9 
Gunn J. E., Peterson B. A., ApJ, 142, 1633 
Haiman, Z., & Holder, G. P. 2003, ApJ, 595, 1 
Hu W., Holder G. P., 2003, PRD, 68, 023001 
Kaiser N., 1992, ApJ, 388, 272 

Kaplinghat M., Chu M., Haiman Z., Holder G. P., Knox L., Sko- 

rdis C, 2003, ApJ, 583, 24 
Kogut A., Spergel D. N., Barnes C, Bennett C. L., Halpern M., 

Hinshaw G., Jarosik N., Limon M., Meyer S. S., Page L. et 

al. 2003, ApJS, 148, 1 61 
Locb A., Zaldarriaga M., 2004, PRL, 92, 211301 
Madau, P., Meiksin, A., & Rees, M. J. 1997, ApJ, 475, 429 
Mo, H. J., Jing, Y. P., White, S. D. M. 1997, MNRAS, 284, 189 
Morales, M. F. & Hewitt, J. 2004, ApJ, 615, 7 
Oh, S. P., & Mack, K. J. 2003, MNRAS, 346, 871 
Oh, S. P. 1999, ApJ, 527, 16 

Pen U., Wu X., Peterson J. B., 2004, astro-ph/0404083 

Santos, M. G., Cooray, A., Haiman, Z., Knox, L., & Ma, C.-P. 

2003, ApJ, 598, 756 
Santo s M. G., Coor ay, A., Knox L., 2004, ApJ submitted, 

astro-ph/0408515 
Scoccimarro R., Zaldarriaga M., Hui L., 1999, ApJ, 527, 1 
Scott, D., & Rees, M. J. 1990, MNRAS, 247, 510 
Shaver P. A., Windhorst R. A., Madau P., de Bruyn A. G. 1999, 

A&A, 345, 380 
Spergel, D. N., et al. 2003, ApJS, 148, 175 
Venkatesan A., Giroux M. L., Shull J. M. 2001, ApJ, 563, 1 
Zaldarriaga, M., 1997, PRD, 55, 1822 

Zaldarriaga, M., Furlanetto, S. R., Hernquist, L. 2004, ApJ, 608, 
622 

Wyithe J. S. B., Loeb A. 2003, ApJ, 588, L69 



REFERENCES 

Barkana, R. & Loeb, A. 2001, Physics Reports, 349, 125 
Bennett C. L., Halpern M., Hinshaw C, Jarosik N., Kogut A., 

Limon M., Meyer S. S., Page L., et al. 2003, ApJS, 583, 1 
Bharadwaj S., Ali S. S. 2004, MNRAS, 352, 142 
Bha radwaj S., Pand ey S. K. 2004, MNRAS submitted 

lastro-ph/ 0410581 1 
Cen R. 2003, ApJ, 591, L5 

Chen X., Cooray A., Yoshida N., Sugiyama N., 2003, MNRAS, 
346, L31 

Chen X., Miralda-Escude J. 2004, ApJ, 602, 1 
Ciardi, B. & Madau, P. 2003, ApJ, 596, 1 
Cooray A. R., Hu W., 2000, ApJ, 534, 533 
Cooray A. R., Hu W., 2001a, ApJ, 548, 7 
Cooray A. R., Hu W., 2001b, ApJ, 554, 56 
Cooray A. 2004a, PRD, 70, 023508 
Cooray A. 2004b, PRD, 70, 063509 
Cooray, A. & Sheth, R. 2002, Phys. Rep., 372, 1 
Cooray A., Furlanetto S. 2004, ApJ, 606, L5 



© 0000 RAS, MNRAS 000, 000-000 



